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Abstract: A analysis is performed to test the resolving power of two-dimensional 
pion interferometry using for illustration the preliminary E802 data on Si + Au at 14.6 
AGeV/c. We find that the resolving power to distinguish two decoupling geometries of 
■ different dynamical models is enhanced by studying the variation of the mean per 

2 . degrees of freedom with respect to the range of the analysis in the qtjQl plane. The 

I preliminary data seem to rule out dynamical models with significant uj,r] resonance 



in 



X 



formation yields. 



^ , In this letter, the results of a new multi-dimensional pion interferometric analysis of 

I preliminary E802 vr^vr" correlation data|l|]-P| on central Si + Au reactions at 14.6 AGeV/c 
are reported. To study the resolving power of multi-dimensional - in contrast to the projected 
one-dimensional interferometry - we consider two dynamical scenarios that predict similar 
correlation functions even though the underlying decoupling geometries differ considerably. 
We use the non-ideal inside-outside cascade model of Refs. . In one scenario long lived 
resonances are neglected while in the other the resonance fractions as taken from the Lund 
model. At higher SPS energies[Q, we showed that long lived resonance production predicted 
by the Lund Model led to similar correlation functions to those expected from models with 
quark-gluon plasma formation @], since long lived resonances could mimic the long slow burn 
of a plasma. Therefore another goal of the present study is to see if a similar ambiguity could 
exist at the lower AGS range. There is some fragmentary data0] indicating that resonance 
production could be suppressed at the AGS, the Lund model predicts similar resonance 
yields at the two energy ranges. We show below that the two-dimensional interferometric 
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analysis has sufficient resolving power to rule out significant long lived resonance production 
in the AGS domain, if the preliminary data are confirmed. 

Recall that under idealized conditions the correlation function, C2(ki, k2), of two identi- 
cal pions probes the decoupling or freeze-out space-time distribution of pions, p{x), through 
C2(ki,k2) = 1 + \p{ki — k2)\'^. However, in actual high energy reactions, final state interac- 
tions, correlations between coordinate and momentum variables, and resonance production 
distort this ideal interference pattern due to Bose-Einstein symmetry (see e.g. [Q-|]TU[). The 
correlation function in this more general case can be expressed as[Q 

where T(g) = {Qc/ q) / {^'^"^'^ is the Gamow factor that distorts the pattern due to Coulomb 
effects, with Qc = 27iam and q = (— (fci — A;2)^)^^^. The complex amplitude, G{ki,k2), 
is the smoothed Fourier transform of the decoupling phase-space distribution D{x,p) = 
{S'^{x — Xf)5^{j) — pf)) , where {xf,pf) are the phase-space coordinates of the pions at their 
decoupling times Xj. For minimum uncertainty Gaussian packets, 

G{k,, k^) = J d^xd^p D(x,p)e'«"e^'^"'/2g(p-i^)V2Ap2 ^ ^^iqx; ^-Kp^/^p^s^ ^ ^2) 

where Ap is the rms momentum width of the packets. The interference pattern depends 
then not only on the break-up phase-space distribution, D{x,p), but in general also on the 
quantal production dynamics characterized here by Ap. The standard notation q = ki — k2 
and K = {ki + k2)/2 is used throughout. The distortions of the interference pattern due to 
long lived resonances modify the amplitude as [Q. 

G{ki,k2) ^ (^f{n' /r) {1 - iqUr/Try^ exp{iqxr - Kpr/Ap^)) , (3) 

r 

where /(vr~/r) is the fraction of the observed 7r~'s arising from the decay of a resonance of 
type r, which freezes-out with final four velocity u!f. The single inclusive pion distribution 
in this notation is -Pi(k) oc G{k,k). The Lund model in the AGS range suggests that 
fiji'/uj) = 0.16, fiji'/K*) = fiji' /{rj + Tj')) = 0.09, while some data seem to favor smaller 
fractions 0. An aim of the present study is to test if multidimensional pion intereferometry 
can discriminate between these possibilities. 

We fit to the correlation data using the following parameterization of the break-up 
distribution g |[ 

,exp{-^-(l9|:)!-(0!-|}.(B-E,)np.) . (4) 

where r = (t^ — z^)2 is the freeze-out proper time, and rj = |log((t + z)/{t — z)),y = 
I log{{E+pz) I {E~pz)) are the space-time and momentum rapidity variables. The correlation 
between these rapidities is estimated Q from the Lund model to be Arj ~ 0.8. 

The goal of the present interferometry analysis is to extract the rms transverse radius, 
Rt, at decoupling and the rms decoupling proper time interval. At. Note that unlike in 



2 



classical cascade models 1 11, 12 , the transverse momentum in the above model is assumed 
to arise entirely from the finite momentum spread of the pion wave-packets through the 
parameter Ap. The y*, Fc, Ap are fixed by fitting the single pion inclusive dataH] with 



Pi (A;) oc G{k, k)oce ^^e-'"^/^m^^/^ (5) 

where T = Ap^/m is the effective transverse mass slope. We search for the optimal values 
of Rt and At by fixing Arj = 0.8 from @, § and = 0.7, y*^ = and Ap^m = 0.17 GeV 
from the inclusive pion distributions Note that we assume implicitly that the chaoticity 
parameter A = 1 in our analysis because, as we have emphasized in [§,0, fits to data treating 
A as a free parameter distort the geometrical scales in a further uncontrolled way. 

For illustration, we recall that for special kinematics with rriiT = = ttlt and Ay = 
yi — 2/2 small, the correlation function has the approximate analytic form when resonance 
production is neglected 

C(.a.)«l + {l + im|A.W(A,V^^^^^))" . (6) 

This shows that the interference pattern as a function of the rapidity difference depends 
not only on the geometrical quantity of interest. At, but also on the width, Ar], of the 
1] — y correlation as well on the width, Y^, of the finite rapidity distribution and on the 
wave-packets size. 

To compare theoretical correlation functions with data projected onto two of the six 
dimensions, we must compute the projected correlation function 

^ f X / rf^fcirf^fc2A(ki, ks) A2{qT, Ql] ki, ks) 

where Pi and P2 are the one and two pion inclusive distributions, and A2 is the experimental 
two pion binning and acceptance function. It is important to note that because of this 
projection. Coulomb corrections do not in general factor out of the integral. All calculation 
were performed using the Monte Carlo importance sampling method adopted in the CERES 
program^, ^. The acceptance function for the E802 experiment was approximated by 

A2(gT,gL;ki,k2) = Ai(ki)Ai(k2)e(20- 101-021) 

Siqi - \hi - h2\)SiqT - |kT, - kT,|) . (8) 

The angles are measured in degrees and the momenta in GeV/c. The single inclusive distri- 
bution cuts are specified by 

Ai(k) = 6(14 < eiat < 28)Qipiat < 2.2 GeV/c)e(y„,„ > 1.5) . (9) 

To assess the statistical significance of the differences between the fits obtained assum- 
ing resonance and non-resonance dynamics, we compute the goodness of fit on a two- 
dimensional grid in the {qTAL) plane, binned with 5qT = SqL = 0.01 GeV/c. In our pre- 
liminary analysis reported in Ref. [|13|, we compared the square of the difference between the 
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theoretical and experimental correlation functions and found an unexpected ridge of high 
variance along the qx = qi diagonal. However, Zajc[jl4| pointed out that this feature could 
have been an artifact of statistical fluctuations in ratios of random variables, and suggested 
that the following variable 



(10) 



should be used instead. In eq. (p!OD, JV-,^ is a normalization factor which is chosen to minimize 
the average and depends on the range in the qr, qi plane under analysis. The indices 
i,j refer to the the corresponding qT,qL bins, in each of which the experimental correlation 
function is given by 



Cexpihj) — A/^ X 



ACexpih j) — Cexpih j) ' 



\ [ A{z,j) I ^ \ B{^,J) I 



■ (11) 



The data for the numerator A{i^j) ± AA{i,j) and denominator B{i,j) ± AB{i,j) were 
obtained from R. Morse with an understanding that the data in this form are preliminary 
and subject to further final analysis. We note that the published one-dimensional correlation 
data in [0 are not in a form suitable for our multi-dimesional analysis. Unfortunately, no 
final multi-dimesional correlation data are available at this time, and we use these data 
only for illustration of the method and to emphasize the need to go beyond traditional 
one-dimensional projections. 

Minimization of the average is performed by exploring the parameter space of Rt and 
At and computing the (x^), averaging over a 30x30 grid in the relative momentum region 
Qt, Ql < 0.3 GeV/c. In the vicinity of the minimum, Rtq, Atq, we determine the parameters 
of the quadratic surface 



{x\Rt, At)) = xLn + o^{Rt - Rto? + KAr - At^)^ . 



(12) 



Recall that the average over N bins is a random variable with unit mean and rms width 

a - 



'2/N if the xihJ) ^ire normal random variables with zero mean and unit rms width. 
For large N, the distribution of the mean per bin is P{x^) oc exp[— (x^ — l)^/2cr^]. For 
the = 900 grid under consideration, a ~ 0.047. 

The most direct measure of the goodness of the fit in this case is rifj = |(Xmm) ^^e 
number of standard deviations from unity of the average x^ per degree of freedom. Inserting 
eq. (|12D into the asymptotic form of the x^ distribution for large A^, the likelihood for the pa- 
rameter Rrp to have a value near the minimum is approximately oc exp[— a^(i?T ■~-Rto)^/2o"^]. 
We therefore estimate the error on the radius to be AR ^ {y2[r(3/4)/r(l/4)]a/a}i/2 ^ 
0.7((T/a)^/^. Similarly, the estimated error on the proper time interval is 0.7{a/P)^^'^. The 
results of this statistical analysis are given in Table 1. 

The lego plots of the Gamow corrected data and the corresponding theoretical corre- 
lation functions with and without resonances, corresponding to the optimized geometrical 
parameters, are shown in Figure 1 for the restricted range 0.005 < qT,(lL < 0.125 GeV/c. 
Part a) shows the 2-D data in the above range. Part (d) shows the corresponding experi- 
mental errors, which vary from AC ~ 0.1 over most of this region to AC ~ 0.35 toward the 
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right most corner at small qr and large q^. In part (b) the best fit (with Rt — 4.6 fm and 
Ar = 3.4 fm/c from Tabic 1 over the wider domain QtjQl < 0.3 GeV/c) assuming no reso- 
nance production is shown. In part (c) the best fit with Lund resonance fractions {Rt = 3.1 
fm and At = 1.6 fm/c again over the wider range) is shown. Finally the x^ihj) fo^ the 
non- resonance scenario is shown in (e), whereas that corresponding to the full resonance fit 
is shown in part (f). Similar large fiuctuations close to the edge of small qr and large qL is 
seen in both cases. However, the overall trend suggests that the no resonance fit is closer to 
the data in the region of small q even though the global is about the same for both fits. 

TABLE 1: 2D-x^ Analysis of Pion Decoupling Geometry 



x'{Rt,At) 


No Resonances 


LUND Resonances 


E802 Data Not Gamow Corrected 




1.125 


1.128 


Rto 


4.1± 1.3 


2.5 ± 1.3 


Ato 


3.3± 1.0 


1.4± 1.0 


a 


0.014 


0.013 




0.021 


0.023 


E802 Data Gamow Corrected 




1.098 


1.104 


Rto 


4.6± 0.9 


3.1± 1.3 


Aro 


3.4 ± 0.7 


1.6± 1.0 


a 


0.027 


0.014 


:] 


0.012 


0.023 



The average per degree of freedom depends, of course, on the range of q under analysis. 
We have tested this by varying the range of the analysis to restricted gr, qh domains, ranging 
from a 1 X 1 grid corresponding to 0.005 < gr, ?l < 0.015 GeV/c, to 2 x 2, etc. as shown 
in Figure 2. However, instead of (x^) we show there the number of standard deviations 
from unity of the averaged P^r degree of freedom. For each n x n grid, N = n"^ is the 
number of degrees of freedom and the standard deviation is expected to be cr = \/2/n. The 
strong dependence of the number of standard deviations from unity as a function of the 
range of the analysis is brought out clearly in Figure 2. We see that the optimal fit including 
Lund resonances is much worse than the fit without long lived resonance distortions when 
the analysis is restricted to the domain qL,qT < Qmax ~ 100 MeV/c, where the correlation 
function deviates significantly from unity. The two models yield similar accepatable fits 
in terms of x^ for Qmax > 100 MeV/c because in that large domain both predict models 
trivially predict nearly unit correlation functions. Also, the mean (x^) begins to deviate 
strongly from unity for both models when data close to the edge of experimental acceptance, 
near qL,qT ~ 0.3 GeV/c, is included. 

We conclude that multi-dimensional analysis has high resolving power in the domain 
of physical interest, and if these data are confirmed, the present analysis would rule out 
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long lived resonance production models in this energy range. However, we note that recent 
preliminary analysis in ref. seems to give systematically lower correlation function values 
in the small g^, qi domain, in which case the resonance or an alternate long lived source 
geometry would be needed. In any case, the results show that 2-D analysis can be used to 
improve considerably the resolving power of interferometry by amplifying in a quantitative 
way the differences between the data and calculations. We note that recent comparisons of 
interferometric data with cascade models, RQMD]TT| and ARC[|T2|, indicate good agreement 
wit one-dimensional projected correlation data. However, the global is quite poor. The 
study of the dependence of that as a function of the range in those cases should help 
resolve where the models fail in a more quantitative fashion. Only in that case can one begin 
to assess whether the data are in fact consistent with standard hadronic transport models or 
require novel collective dynamics. In general, the analysis could be sharpened considerably 
by independent measurements of long lived resonance abundances to reduce that source of 
interferometric distortions. 
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Figures Caption 

Fig.l Negative pion correlations in central Si+Au reactions []I[] as a function of transverse 
and longitudinal momentum difference Qt-, Ql- The preliminary E802 data corrected 
for acceptance and Coulomb effects are shown in part (a). Parts (b) and (c) show 
theoretical correlation functions filtered with the E802 acceptance for cases without and 
with resonance production, respectively. The transverse radius and proper time interval 
characterizing the pion decoupling distribution were obtained by minimizing the over 
a wider region with gr, Q'l < 0.3 GeV/c. Part (d) shows the distribution of X^{<1t,(1l) 
for the case with resonance production in part (c). The average parameters are 
listed in Table 1. 



Fig. 2 The number of standard deviations from unity of the P^r degree of freedom as a 
function of the range Qmax of the analysis in the qt, Ql plane. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9501035vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9501035vl 
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